ON PLOUFFE'S RAMANUJAN IDENTITIES 



LINAS VEPSTAS 



ABSTRACT. Recently, Simon Plouffe has discovered a number of identities for the Rie- 
mann zeta function at odd integer values. These identities are obtained numerically and 
are inspired by a prototypical series for Apery's constant given by Ramanujan: 

77T 3 ~ 1 



C(3) = ^--2V 



180 " ^ n 3 (e 2 ™ - 1) 

n— 1 v ' 

Such sums follow from a general relation given by Ramanujan, which is rediscovered and 
proved here using complex analytic techniques. The general relation is used to derive 
many of Plouffe's identities as corollaries. The resemblance of the general relation to the 
structure of theta functions and modular forms is briefly sketched. 



1. Introduction 

Inspired by an identity for £(3) given in Ramanujan's notebooks [4, chapter 14, formu- 
las 25.1 and 25.3], 



7-7T 3 °° 1 

« 3 > = ^- 2 £ 



180 " ^ n 3 (e 2 ™ - 1) 

n— 1 v ' 

Plouffe describes a set of similar identities] 8 1, [9] that were discovered numerically using 
arbitrary-precision software. For example, Plouffe gives an identity for ((7): 



v ' 56700 ^ n 7 (e 2ml - 1 



(e 2 ™ - 1) 

n— 1 v ' 

This text provides an analytically derived formula for expressions of this type. The result- 
ing general formula, valid for integer m > 1, is 

C(4m-1) = -2^ 



i( 4m-l ( e 27rn _ ]_) 
n 2m 



B2j Bi m -2j 



(2j)!(4m-2j)! 

where Bk is the /c'th Bernoulli number. 

Plouffe also notes similar relations for the other odd integers; for example, 

_ 7T 5 72 " 1 2 " l 



294 35 ^ n 5 (e 2 ™ - 1) 35 ^ n 5 (e 2 ™ + 1) 

n—1 v / n—1 v 7 
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The general form for this type of expression may be shown to be 



[l + (-4) m - 2 4m+1 ] C(4m + 1) 




4m+l ( p 2iTn 



+2 [2 4m+1 - (-4)™] 

n=l 





The methods described in this text also allow for a large generalization of these types of 
sums. Defining 



these generalizations follow from a modular equation relating (r) to (— 1/ r) for odd 
integers k, the derivation and proof of which is the one of the main topics of this note. The 
modular relation is not new; it appears in Ramanujan's Notebooks [4, Chapter 14 Entry 
21] as 



where a > 0, (3 > with a/3 — tt 2 and n any positive integer. Berndt implies that this 
formula is the most studied of all the notebooks; it has been independently discovered per- 
haps a half-dozen times, and proven twice as often. It has been generalized to L-functions, 
and to rational values of k; Berndt provides a long list[4 | of the various proofs and gener- 
alizations made. 

Much of this paper is devoted to (yet another! independently discovered) proof of this 
relation, followed by a series of lemmas that provide the connection to Plouffe's results. 
In searching for curious and interesting special cases of this relation, one senses that only 
the tip of the iceberg has been seen. Unexplored possibilities include, for example, consid- 
ering t S Q[i], the field of Gaussian rationals, or from considering Diophantine roots of 
quadratics. 

The rest of this paper is roughly laid out as follows: The second section provides a re- 
view of previous related results. The third section gives a relationship between the sums 
and the polylogarithm, and thence to an integral on the complex plane. The fourth sec- 
tion examines the related contour integral, which is easily integrated via Cauchy's residue 
theorem to give a finite sum involving the Bernoulli numbers. The fifth section relates the 
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contour integral to the polylogarithm integral, thus resulting in a functional equation for 
Pfc(r). The sixth section applies the functional equation, providing a variety of lemmas, 
many of which explain Plouffe's discoveries. The seventh section give a pair of relation- 
ships on the Bernoulli numbers that arise naturally in this context. The eighth section 
explores the modular nature of the relations on -Pfc(r), followed by a conclusion. An ap- 
pendix gives a derivation of an integral representation of the polylogarithm, that is central 
to the analysis. 

2. Related sums 

A large number of similar sums have been explored before; this section reviews some 
of these. Perhaps the most forthright is a sum given by Ramanujan in a famous letter to 
Hardy|6|, stating that 

l 13 2 13 | J*_ 1_ 
A generalization of this sum, 



E 



„4fc+l R 

n -c>4fc+2 



'-1 4(2fc + l) 

is proved by Berndt[2|, and attributed to Glaisher[5]. This and many related results are 
derived by Zucker flT2l . based on the theory of Jacobian elliptic functions. A similar result 
is stated by Apostol in the form of an exercise[l | see exercise 15 at end of chapter 1.]: 

°°. „4fe+l o4fc+l -I 

TT^ ~ 8fc + 4 Bik+2 

n=l-n Odd 

Many sums resembling those in this note are given by Zucker lfTTIl . Some of these are 

1 coth(7rmx) — 1 



m—l 



rn 



n fornix _ J) 2 ^ 
n— 1 m—l 



m 



It 



1 — tanh(7rma;) 



rn 



^n(e 2 ™ x + l) 2 

ri— 1 v y m—l 

(-!)" _ 1 f 1 - tanh(7rroa;) 



^ n(e 2 ™ x + 1) 2 ^ v 

n— 1 m—l 

The fount of inspiration for such sums is Ramanujan. Sums given in Chapter 14 of Part II 
of the Ramanujan's Notebooks[4| include entry 8: 

sinh(2cmfc) sa\(2Bnk) a , . ,, 3 ,„,„ k 

n— 1 n— 1 

and another similar one relating cos and cosh. Above, fc is any positive integer, ad — ir 
and < 3k < tt. 

Entry 1 3 generalizes the sums mentioned previously, 

« k E - (-0* E = - (-m ^ 
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This time, one takes a/3 = 7r 2 and k > 1 an integer. The above follows from sums on the 
divisor function, as is frequently noted. 

Sums involving pairs of Bernoulli numbers also appear in the analysis of the Dedekind 
eta function. Thus, Sho Iseki's transformation formula, as described by ApostolQ] see 
theorem 3.5], is 

wo\ i /, a l\ ^ ( 2 \ B 2 _ n {a)B n {(3) 
A(a, p, z) = A I 1 — p, a, - I — ttz > 

where A is given by 



n J (iz) n 



A(a, 0, z) = E [K z i r + a ) ~ ifi) + K z i r + 1 - a) + iP)] 

r=0 

and 



m = E 



m— 1 

A sum linking the Bernoulli and Euler numbers is given by BerndtQ: 



" l S ech[(2n-l)7rV3/2 
^ ( ' (2n-l) 6k + 1 

n=l v ' 



Y /( f 2m +\, „f 6fc - 2m „ cos [(2m + 1) 
2 y ' ^ (2m + 1)! (6k- 2m)! L v ; 



Perhaps the results that are closest to those presented in this paper are those noted by 
Borwein, et al flT] Section 5], and in particular, given a very similar result involving ((Am— 
1) and C(4m + 1). 



3. The Polylogarithm 
The recurring theme in Plouffe's identities is the sum 



S s( x ) - zl n s (gin _ 



with s usually an odd positive integer and x = ir or x — 2ir or possibly other interesting 
values, such as x = n^/m for some integer m. This sum may be converted into a sum over 
polylogarithms, and subsequently into an integral. The integral may, after some difficulties, 
be converted into a contour integral, whereupon it may be evaluated by Cauchy's residue 
theorem. The result is a finite sum whose general structure resembles those of Plouffe's 
and Ramanujan's identities. This section develops the first part of this analysis. 
To find the polylogarithm, one expands 



= EE 

n— 1 m— C 

oo 

- E Lis 



x XJ e -xn(m+l) 
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which is generally valid for ^Rx > 0. The above is easily obtained by applying the expan- 
sion 



oo 

z r 



l-z ^ 

m=0 

and the definition of the polylogarithm: 



71=1 

The polylogarithm may be expressed in terms of an integral as 

pc-\-ioo 



i rC-\-lOQ 

Li s (e-") = — / r(z)((z + s)u- 



z dz 



(See Appendix for derivation). Here, T(z) is the Gamma function. The line of integration 
is taken to be to the right of all of the poles in the integrand, namely, c > 1. Using this in 
the summation, one obtains 



1 pc 

S *( x ) = I r(z)((z + S )(xmy z dz 



1 rc+ioc p/ x 

«L -^«' + •>«'>* 



The exchange of the order of summation and integration is justified precisely when one has 
c > 1. The last integral has poles at z = 1 and z + s = 1 coming from the zeta functions 
and poles at all of the non-positive integers coming from the Gamma function. 

If the integral can somehow be converted into a closed contour on the left, then it may 
be evaluated in a straight-forward way by means of Cauchy's residue theorem. Performing 
this closure is in fact harder than one might hope, as there are non-zero contributions to the 
contour from its closure. The next section evaluates the Cauchy integral, assuming that the 
contour can be closed. The section after that computes the contributions from closing the 
contour integral. 



4. The Contour Integral 
Define the contour integral as 

I s (x) = ±J^t(z + s)«z)dz 

where the contour 7 encircles the poles atz = l, z + s = 1 and z = 0, — 1, —2, ... in the 
usual, right-handed fashion. Then, one uses Cauchy's theorem, which states that 

2-ki J z — a 

for simple poles, and that 

''<•'- 55/ (TV 

for double poles. For the pole at z = 1, one obtains the residue 

Res(z = 1) = «i±l) 
x 
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For the poles at z = — n, one obtains the residue 

(- x y 



Res(z = —n) = - — -j—((s — n)((—n) 



and so one has 

I s (x) — Res(z = 1 — s) + 
For s not an integer, one has 



C(* + i) 



+ E 



n=0 



(—x) n 

i-J-C(«-n)C(-n) 



Res(z = l-s) = ^T^C(1 - s) 

However, the interesting case is for s = k a positive integer. In this case, the pole overlays 
another pole from the Gamma, and one has a double pole. This is just a little trickier to 
evaluate: 



1 



z=l-k 



d_ 

dz 



-dz 



--l-k (z + k-1) 
(z + k-l) 2 T(z) C{z + k) 



CM 



z=l-k 



To perform the derivative, one will need to use the identities 

d 



ds 



(s-l)C(s) 



= 7 



8 = 1 



where 7 = 0.577 ... is the Euler-Mascheroni constant, and 



r(n+l) 

where ip(z) is the digamma function, and H n is the n'th harmonic number. Putting these 
together, one obtains 



2m 



^t(z + s)t(z)dz = [C'(l -k) + (H fc _! - In 2n) ((1 k)] 



Adding this to the other contributions, one gets 

C(k + 1) 



h{x) 



+ 



E 

n = 

n 7^ k — 1 



(—r) n 

i-J-C(fc-n)C(-n) 



+ 7fc^i)i" [c ' (1 _ fc) + ~ ln27T) c(1 _ fc)] 

When k is an odd integer, the above simplifies in two ways. First, ((1 — k) vanishes, 
because the zeta function vanishes at all negative even integers. Similarly, the infinite sum 
becomes finite: when k is an odd integer, one has either Q{k — n) — or ((— n) = for 
all n > k. Thus, for k an odd integer, one has 

= £ t^C(*-„)C(-n) + fc^C(l-*) 



n = 
n ^ fc - 1 



(fc-1)! 
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For the remainder of the paper, it is assumed that, in this context, k > 1 is an odd integer, 
unless explicitly stated otherwise. The above evaluation of the contour integral is the main 
result of this section. To see that this is a key result, one may substitute k = 3 and x = 2n 
to obtain 

which should be recognizable as a portion of Ramanujan's identity. For k — 7, one has 

I9ir 7 



/ 7 (2tt) = -C(7) 



56700 

which resembles one of the results given by Plouffe. To complete the connection, one must 
relate the contour integral Ik (x) to the sum S k (x). This is done in the next section. 

First, however, to drive the point home, one must observe that most of the terms in the 
above expression are rational multiples of powers of it. This follows from the zeta function 
being related to the Bernoulli numbers B n at even integers: 

_ 1 (2n)^B 2n 



2(2n)\ 

for integer n > 0. At the negative values, one has 

B n +i 



C(2n) = (-1)' 
values, one 
C(-«) = -- 



n+l 
while the derivative is 

C'(-2n) = (-irJ^C(2n + l) 

for integer n > 1. Using these in the above expression for Ik{x), and rearranging terms a 
bit, one gets 

fe-il 



-2I k (x) = C(k) 



i-f— V 



The above introduces the binomial coefficient 

n \ n\ 



k J k\{n-k)\ 

Although the imaginary number i = \f~^\ appears in the above, it is always squared, and 
thus is just a sign-keeping device. Every term in the sum is purely real. 

When the contour integral is written in this form, it may now be seen that for x being 
any rational multiple of n, that is, x = pn/q for any integers p, q, that the coefficient of 
£(fc) is a rational number, and that the second term is another rational times ir k . 

A further curiosity in this regard is noted by Plouffe: if one takes x — ^\fp~fq for 
integers p and q, one also gets simple expressions: because k is odd, the coefficient of 
is still a rational, and the coefficient of n k is \Jpjq times some rational. For rational x, one 
still has that the sum is a rational polynomial in ir 2 , and for x — \fwp~Jq, one still has that 
the sum is a rational polynomial in n. The ocean-full of rationals here suggest that some 
sort of p-adic analysis might be interesting. The appearance of the square root suggests 
that there is a relation to complex multiplication, or that one may have interesting results 
on the field of Gaussian integers. 
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5. Evaluating the Contour Integral 
The goal of this section is to relate the sum Sk(x) to the contour integral Ik(x). 
Theorem 5.1. For odd integers k, one has that 

5 fe (,)=4(,) + (-l)(- 1 )/ 2 (^) fc - 1 ^(^) 

For the remainder of this text, this will be referred to as the "functional equation for Sk ". 

Proof. To prove this result, consider evaluating the contour integral I s (x) for a tall rect- 
angular contour surrounding the poles at z = 1, 0,-1, . . . , 1 — s. Thus, write I s (x) = 
A + B + C + D with A being the integral from c — ih to c + ih for a constant c > 1 and 
the height h large, eventually taking the limit h —> oo. That is, A forms the right hand side 
of the rectangular contour. In the limit of h — > oo, one has by definition 

A = S s (x) 

Let B and C be the top and bottom of the contour, so that for B, the integral runs from 
c + ih to ih — s — e leftwards. For C, the integral runs rightwards from — ih — s — e to 
c — ih; here we take e > 0. The integral D on the left hand side of the rectangle closes the 
contour, running downwards, from ih — s — e to — ih —s — e. 

The integrals B and C will vanish in the limit of the height h — > oo. This can be easily 
seen after a simple change of variable: 

-i Ah— s—e 

B - 55 L 4^ + ^ 

1 /- s_e T(u + ih) , . , 

-y— C(s + u + ih)((u + ih) du 



2ni 

Much of the integrand is 0(1) in h, or polynomially thereabouts. The integrand is domi- 
nated by the Gamma function, which, from Stirling's approximation, may be seen to be 



T{u + ih) = O (e-^/ 2 ) 



The complex conjugate argument applies to C, and thus B and C vanish in the limit of 
h — ► oo. 

To evaluate D, begin by writing 

d = — / A2 c{s + z)c{z)dz 



2ni 



ih — s — e 
ih 



1 f ln 

= -7—x s+t x u r(-u-s-e)((-u-e)({~u-s-e)du 
2™ J-ih 

after a change of variable z = — m — s — e. One then applies the functional equations for 
Gamma: 

T(l-z)T(z) 



and for zeta: 
to obtain 



sin ttz 

C( s ) = 2 s ^- 1 sin^r(i- s )C(i- s ) 



2¥^V27r; V27r; J- 4 ft ^ 47I . 2y / sin(u + s + e)?r 

T(l + u + e)C(l + u + e)C(l + s + u + e) du 
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Another change of variable, this time asw = 1 + u + etake the integral to a slightly more 
recognizable form: 



D - -— (— / ( x 

m\2n) ./,., ,, V l- 2 



1 ! ' ' x ycos^ cos^tt 



! / sin(w + s)ir 



T(w)((w)((w + s) dw 



Next, by taking s = k to be an odd integer, the trigonometric piece simplifies and looses 
its w dependence: 

COS |7T COS ^7T = l ( _ 1)(fc+ l)/ 2 

sin(w + k)ir 2 

Pulling out this piece, one regains a recognizable integral, so that, in the limit h — > oo, one 
finally reaches the claimed result: 



D 

that is, 



for odd integer k. □ 

The next section will review the application of this and the preceding sections to spe- 
cific, simple values of x, thus regaining many of Plouffe's sums. 

6. Lemmas and Applications 
The first corollary demonstrates Plouffe's simplest sums for ((4m — 1). 
Corollary 6.1. For m integer, one has 

^4m-i(27r) = 2S , 4m _i(27r). 
Proof. Substitute x = 2tt in the functional equation. □ 
This corollary provides the first concrete result of this exposition, namely that 

oo -. ,r. \4 m _i 2m / \ 

% „,„-, (e L_ 1) = -«*» - 1) + g ( t; ) <-•>' ^»-« 

which completely resolves one set of relationships given by Plouffe. The functional equa- 
tion opens additional possibilities. By substituting x = 2np/q, one obtains, for k = 
4m - 1, that 

Thus, for example, by choosing p = 2, g = 1 and k = 3, on obtains 



29tt 3 2 ^ i 

rt=l 



4 1 

+ 



and similarly, for k = 7, 

409tt 7 2 



f . = 409tt^ _2y^2. 
^ ' ~ 94500 5 ^ n 7 

and one may proceed in a similar manner. 



e™ - 1 e 47r " - 1 



+ 
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The next corollary shows that something more is needed for C(5), C(9)> an d so on > since 
the most direct approach does not give any information for such sums. 



Corollary 6.2. For k an odd integer, one has 



Proof. This may be proved by applying the functional equation twice in a row. That is, it 
may be proved by substituting x — > An 2 /x in the functional equation and then employing 
the result. □ 

Corollary 6.3. For m a positive integer, one has 

I 4m+1 (2ir) = 0. 

Proof. Substitute x = 2n in the preceding corollary. □ 

Following directly from the above is another corollary: 
Corollary 6.4. For m a positive integer, one has 

2m+l p 

0= V ( iv" Bim - 2 J+ 2 ^ 

^ > (4m-2j + 2)!(2j)! 
Proof. Write out the value of Ii m+ i (2tt) in detail. □ 

From the above, it should be clear that the functional equation for Sk does not provide 
any statements about Sk when k = 4m +1. To obtain results on sums involving k = 
Am + 1, one must introduce 

n=l y ' 

Theorem 6.5. One has 

T s {x) = S s (x)-2S s {2x) 

Proof. This may be proved by re-writing in terms of the polylogarithm, along the lines of 
the earlier development: 

oo 

T,{x) = -Y t {-\) m U a (e' xm ) 
The even and odd terms are regrouped, as 

oo 

T„{x)= ^Li s (e- xm )-2U s (e- 



m— 1 

which is seen to be a sum of S s 's. □ 



The results for ((5), etc. follow from a critical observation: that 

S s (x + 2ni) = S s (x) 
is a periodic function. This periodicity is employed directly in the next theorem. 
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Theorem 6.6. For positive integer m, one has 

Sim+lfa) = 7 4to+ i(27t(1 + i)) 

74m+l(27T) 



+(-!)' 



2.4 m S 4m+ i(2^)-4 m /4 m+ i(7r) 



4m 

Proof. Using periodicity, one writes, for k = Am + 1, 

S fe (27r + 2m) = S k (2n) = I k (2n(l + i)) + 2^' 2 e^ k - 1 ^ i S k (n(l - i)) 

The series 5fe(7r(l — i)) doesn't have an imaginary part; rather, it is an alternating series, 
which may be expanded and written as 

b k {ir{l - i)) = -T k {ir) H ^ 

The Tk (it) term may be eliminated by writing 

T k (w) = S k (7r) - 2S k (2n) 

and the S k (tt) term may be eliminated by 

SfcW = I k (*)+2 1 - k S k (4n) 

= I k (TT) + 2- k [S k (2TT)-T k (2n)} 

Performing the various substitutions suggested above proves the theorem. □ 

As an example of the application of the above theorem, take m = 1, that is, k = 5. One 
easily finds that 

15 ^ 5 
W = --C(5) 



32^ ' 9-64 
and that 

J 5 (27r(l + i)) = -|c(5) + g^g 

Combining these, one gets 

C(5) = ^-^[T 5 (27r)+36S 5 (2^)] 

which is given by Plouffe. The theorem may be used to generate similar expressions for 
allC(4m+l). 

Curiously, the theorem yields results for m = as well. In this case, one finds 

Si (270 +11(270 = J-7log2 
6 4 

Many identities for n are possible by taking two different expressions for a given zeta, 
and subtracting them, leaving behind a rational combination of the sums and n. Thus, for 
example, the following theorem for Apery's constant: 

Theorem 6.7. A series expression for n 3 is given by 

7T 3 = 720 • S 3 (w) - 900 • 5 3 (2tt) + 180 • 6*3(477) 
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Proof. This follows by taking the general expression for k = 4m — 1: 

16 ro 5 4m -i(7r)+45 4m -i(47r) = 16 m / 4m -i(7r) 

= -^C(4m-l)[16 m + 4] 



2m R R 

-(2ir) 4m - 1 Vf-4V— ^ 
( } ^ ( } (4m-2j)!(2j)! 



solving for ((4m — 1) and then using 



2m 



2S 4m -i(2^) = / 4m _i(27r) = -C(4m - 1) + (2^) 4 ™- 1 ]T (-l) j 
to eliminate the appearance of the zeta. The resulting expression may then be solved for 



(4m -2j)\ (2j)\ 



Relationships involving square roots also arise naturally. 
Theorem 6.8. One has 

5tt 3 

~72~ 2 



□ 



C(3) 



\ S 3 faVs) 



Proof. This follows from the general expression given earlier, that 

fc-i 



Sk 



(2-K-p 
V 9 



Sk 



f2irq 
V P 



Here, making the substitution q = ^/p one obtains 



5 fc (2^)+P (fe - 1)/2 ^(^) 



-C(fc) l - (-p)^- 1 )/ 2 



(fe+l)/2 R 



2j 



(fc + l-2j)!(2j)! 



The specific result follows after choosing k = 3 and p = 3. 

7. Some Bernoulli number identities 
In addition to the previously noted identity 



□ 



2m+l 



0= £(-!) 



-B4m-2j+2 



2j 



J=0 



(4m-2j + 2)! (2j)! 



there are several other identities on sums of Bernoulli numbers that result from the previous 
developments. These are briefly stated here. 



Theorem 7.1. For integer m, one has 



= £(-4)' 

3=0 



B 



4 ^r + 2- 



B. 



Am—Aj 



B 



Aj+2 



_ (4m - 4j + 2) ! (4j) ! (4m - 4j) ! (4j + 2) ! _ 
Proof. Consider the sums resulting from = I 4m+ i(27r(l + i)) — J 4m+ i(27r(l — i)). □ 
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Theorem 7.2. For integer m, one has 



2m 



E, i\fc B 4m-2k B 2 k .sj B*m-ij Bjj 

( =o { ' (4m-2fc)! (2k)\ " " ^ ' (4m -4j)! (4j)! 

m-l p 

; (4m-4j-2)!(4j + 2)! 



'4m-4j-2 £>4j+2 

3=0 

Proof. Consider the sums resulting from the identity 

/ 4 m-i(27r) = / 4ro _i(27r(l + »)) + / 4m _i(27r(l - i)) 

□ 



8. Modular Relations 

The various sums and quantities above can be seen to be quasi-modular by making 
a simple change of variable, namely by making the substitution x = 2nir. By "quasi- 
modular", it is meant that the various terms almost have simple behaviors under the Mobius 
transformation r — ► (ar+b) / (cr+d) for integer a, b, c and d. In this respect, the sums bear 
close resemblance to theta functions, which obey similar relations. These relationships are 
brought to focus here. 

First, define K k (r) — I k (2iriT). This change of variable results in a definition which 
seems simpler than that for I k . An expansion in r is often referred to as a "Fourier series" 
in the context of hyperbolic geometry: 

K k (r) Tfc ' 1 - 1 C(^ - ^ ( V /2 r 2 ^ Bk+1 ^ 
Kk[T) ~ 2 C(fc) 2r 2- T (2j)!(fc + l-2j)! 

This quantity is almost a modular form of weightfc — 1, in that 

K k (—) = -T*- k K k { T ) 



T 

Its only "almost" a modular form, because it is not periodic in r, that is 

K k (r + 1) ? K k (r) 
By contrast, Pfe(r) = S k (2iriT) is periodic: 

P fe (r + 1) = P k (r) 
but is not quite modular under inversion: 

One may define a simple variant of the sums that does have a simple transformation under 
inversion, namely 

M k {T)=P k ( T )- l -K k (T) 

which transforms as 



T 



However, M k is then not periodic. 
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Figure 8.1. Phase of P_ i (g) on the unit disk 




This graphic shows the phase 



argP_i( g ) = arg^ — — - 

71 — 1 

where arg f(z) = 3 log f(z) is the usual arg of a function. The color scheme is such that 
black represents areas where arg > and red represents areas where arg <0. The absence 

of other colors indicates that the phase is rather closely confined to the vicinity of for 
most all of the disk. Numerically, the absolute value of the phase is smaller than 10~ 3 for 
much of the disk. In particular, this indicates that there are no zeros at all in the interior of 
the disk, as a zero would be surrounded by a region where the phase wraps around by 2n. 
The function P s (q) does have poles at q — e 27r,m /' i for all rationals m/n; these are 
visible at the edges of the disk. The fractal nature of this image is the characteristic 
signature of a modular form of weight 2; the self-similar regions are just copies of the 
fundamental region of the modular group SL(2, Z). 
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Figure 8.2. Phase of P-s(q) on the unit disk 




This graphic shows the phase 



oo r 

x n 
arg P_ 5 (<z) =arg2^— — - 

71 — 1 

The color scheme is as in the previous image, and similar conclusions apply: there are no 
zeros at all in the interior of the disk. The absolute value of the phase is tiny: numerically, 

it is within 1CP 9 of zero for much of the disk. 
The fractal nature of this image is the characteristic signature of a modular form of weight 
6; it should be compared to the image of the modular invariant g% shown in the next 

image. 



The analytic structure of Pk is curious: it has a pole at r = ioo and. by the inversion 
formula and periodicity, at every rational value of r. This is clearly visible in the graphic 
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This figure shows the imaginary part of the modular invariant (73, one of the invariants of 
an elliptic curve. Specifically, it shows the imaginary part of 



8^ 6 

9a{q) = -27- 



1 - 504 



\ " 



n 5 q" 



which is a modular form of weight 6. The colors are chosen such that black represents 
areas that are negative, blue and green represent areas with smaller values, and red those 

areas with the largest values. 
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below, which shows P3011 the (/-series or "punctured disk" coordinates q = e 27TlT : 

Psiq) = n s ( g » _ 1) 

n— 1 x 7 

In this form, the relation to modular forms becomes the most immediate. Consider the 
Lambert series: 

00 1 n 00 

£ = " c(s) " Ps{q) = S a - s(m)<r 

n— 1 m— 1 

Here, <r s (m) is the divisor function: 

tr s (m) = n s 

with the notation n\m denoting that the sum extends over all divisors n of m. This should 
be compared to the Eisenstein series [ 1 see section 3.10] 

G 2k (q) = 2C(2fc) + J^-l- ^k-i(m)q m 

^ ' ra—l 

etc. In essence, we are dealing with an extension of the Eisenstein series to general complex 
argument. XXX unfinished; this should be elaborated on a bit more. XXXX 

9. Conclusions 

Given any particular rational and set of sums, there is a rich variety of possible rela- 
tionships that can be formed. It would be interesting to characterize the number of these 
relationships as a function of the number of sums involved. Thus, for example, given only 
63(271-), there seems to be only one relationship that can be formed to give C(3). However, 
given Ss(tt), 5*3(271) and 5*3 (47r), there seem to be two distinct relationships, of which a 
linear combination can be taken: thus there is a one-dimensional linear space of solutions. 
It would be interesting to know what the higher order combinatorial possibilities may be, 
the dimensionality of the spaces they form, and how these depend on the selected series. 

Another interesting avenue of research would be to strengthen the connections to mod- 
ular forms, to see if additional insights can be gained by following that route. A different 
and fruitful avenue for understanding almost-modular structures is to pose them in terms of 
functions on a one-dimensional lattice (in the sense of the Ising or Potts model of physics), 
where translation is described by the shift operator. It is known that such lattice models are 
isomorphic to measure-preserving dynamical systems, such as the Baker's map or more 
generally subshifts of finite type. In particular, it seems possible to re-express the func- 
tional equation for Sk, as given in section 4, in terms of (logarithms of) eigenfunctions of 
the Baker's map. 

Thanks to Simon Plouffe for generating interest in such sums. 

10. Appendix: Polylog Integral 
This appendix proves the following theorem: 
Theorem 10.1. The polylogarithm may be written as the integral 




T{r)Q{r + s)u r dr 
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Proof. The proof below is cribbed from the Wikipedia article on polylogarithms| 10|. One 
begins by writing the the Mellin transform of the polylog, as 



M s (r)= r Us (ye- u )u r - 
Jo u 



Using an integral representation of the polylogarithm 

i r°° t 
r ( s ) Jo w ~ le 



Li * H = — / — Tzr^r dt 



and substituting, one obtains 

MJr) = —- / / — — dtdu 

A change of variable t = ab and u = a(l — b) with dtdu = a da db gives 

i fl roo „r+s—l 

M s (r) = — V-\l-bY- 1 dbl — a da 
r(s) Jo Jo V 1 e Q - 1 

= r(r)Li r+s (y) 
The inverse Mellin transform may now be employed to write 

i />c-Hoc 

Li s (ye-") = ^-. u- r T{r) Li r+S (y) dr 

By setting y = 1, one then uses Li s+r (l) = £(s + r) to obtain the desired result. □ 
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